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Abstract. We present a simple, but efficient, way to calculate connection 
matrices between sets of independent local solutions, defined at two neighboring 
singular points, of Fuchsian differential equations of quite large orders, such 
as those found for the third and fourth contribution (x' 3 ' and x' 4 ^) to the 
magnetic susceptibility of the square lattice Ising model. We deduce all the 
critical behavior of the solutions x' 3 ' and x^ A \ as wen as the asymptotic behavior 
of the coefficients in the corresponding series expansions. We confirm that the 
newly found quadratic singularities of the Fuchsian ODE associated with 

x (3) 

are not singularities of the particular solution x^ itself. We use the previous 
connection matrices to get the exact expressions of all the monodromy matrices of 
the Fuchsian differential equation for x' 3 ' (and x' 4 ' ) expressed in the same basis 
of solutions. These monodromy matrices are the generators of the differential 
Galois group of the Fuchsian differential equations for x' 3 ' (and x' 4 '); whose 
analysis is just sketched here. As far as the physics implications of the solutions 
are concerned, we find challenging qualitative differences when comparing the 
corrections to scaling for the full susceptibillity x at high temperature (resp. low 
temperature) and the first two terms x^ 1 ' and x' 3 ' (resp. \^ and x' 4 ^) • 
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1. Introduction 

Since a pioneering, and quite monumental, paper on the two-dimensional Ising 
models, it has been known that the magnetic susceptibility of square lattice Ising 
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model, can be written as an infinite sum of (n — l)-dimensional integrals EJ 0J 
0013 contribut ions: 

oo 

X(T) = Ex W m (1) 

n=l 

The odd (respectively even) n correspond to the high (respectively low) temperature 
domain. These (n— l)-dimensional integrals are known to be holonomic, since they are 
integrals of holonomic (actually algebraic) integrands. Besides the known x and 
%( 2 ) terms, which can be expressed in terms of simple algebraic or hypergeometric 
functions, it is only recently that the Fuchsian differential equations satisfied by 
the an d X terms have been found |SJ UJ E3- These two exact differential 

equations of quite large orders (seven and ten) can be used to find answers to a 
set of problems traditionally known to be subtle, and difficult, for functions with 
confluent singularities, like the fine-tuning of the singular behaviors for all the 
singularities (dominant singular behavior, sub-dominant, etc.), accurate calculations 
of the asymptotic behavior of the coefficients, etc. 

Recall that the third, and fourth, contribution to the magnetic susceptibility x^ 3 \ 
and x^ 4 \ are given by multi-integrals and each is, thus, a particular solution of the 
corresponding differential equation. These differential equations exhibit a finite set of 
regular singular points that may (or may not) appear in the physical solutions x 
and X ■ Besides the physical singularities and the non physical singularities s = ±i 
(where s = smh(2iT), K being the usual Ising model coupling constant, K — (3 J), 
it is commonly believed that the X s have, at least, other non physical singularities 
given by B. Nickel jHlCj- The dominant singular behaviors at all these (non physical) 
singularities (x' 3 '' and x' 4 -*) have also been given by B. Nickel. The differential 
equations of the X^'s, which "encode" all the information on the solutions and their 
singular behavior, in fact, allow us to obtain not only the dominant, but also all the 
subdominant singular behavior, hardly detectable from straight series analysis. It is 
thus of interest to get (or confirm) these singular behaviors from the exact Fuchsian 
differential equations that we have actually obtained for x^ and X^ and, especially, 
the singular behavior at the two new quadratic singularities, 1 + 3w + 4w 2 = 0, 
(where w = s/(l + s 2 )/2) found for x (3) |H|- 

The physical solution x^ 1S defined by a double integral on two angles and is 
known as a series obtained by expansion (then integration) of the double integral at 
w = (or s = 0). It is certainly not simple to obtain the x^ expansion around 
(say) the ferromagnetic critical point w — 1/4, due to a singular logarithmic behavior. 
However, one can overcome this difficulty since, with a differential equation, it is 
straightforward to obtain the formal series solutions at each regular singular point 
(i.e., a local basis of series solutions). By connecting the formal solutions around 
w = and the formal series solutions around another regular singular point like 
w = 1/4, one will be able to express the particular solution x^ (and also all the other 
formal solutions) as a linear combination of solutions valid at w = 1/4. The seven local 
solutions at w = will, then, be given by the product of a 7 x 7 matrix with the vector 
having the seven local solutions at w = 1/4 as entries. In other words, succeeding in 
obtaining these connection matrices amounts to building a common (global) basis of 
solutions valid for all the regular singular points. Furthermore, with these connection 
matrices, we obtain, in fact, the analytic continuation in the whole complex plane of 
the variable w, of x^ and x^\ which are known as integral representations. 

Note that, remarkably, the Fuchsian differential equation for x^ has simple 



Galois group for the Ising model 



3 



rational, and algebraic, solutions. These rational or algebraic solutions, known in 
closed form, can be understood globally. One can easily expand such globally defined 
solutions around any singular point of the ODE, and follow these solutions through 
any "jump" from one regular singularity to another one, and, therefore, from one well- 
suited basis to another well-suited basis. For a function not known in closed form, 
like the "physical" solution the decomposition on each well-suited local basis 

associated with every singular point of the ODE, is far from clear. The correspondence 
between these various (well-suited) local bases associated with each singular point 
of the ODE, is typically a global problem and, thus, a quite difficult one. One 
clearly needs to build effective methods to find such connection matrices in the case of 
Fuchsian differential equations of order seven, or ten (x*- 3 ) and X )i or of much higher 
orders (x^, X > etc.). With a method of matching of series, we will show that the 
connection matrices matching these various well-suited bases of series-solutions can 
be obtained explicitly. The entries of these matrices can be calculated with as many 
digits as we want. We will show that we can actually find the exact expressions of 
these entries as simple algebraic expressions of (in the case of the Fuchsian ODE's of 
X^ 3 ' and x'" 4 '') powers of tt, hi(2), ln(3) and various algebraic numbers or integers, 
together with more "transcendental" numbers like the "ferromagnetic constant" 1$ 
introduced in equation (7.12) of [Q: 

1 f°° f°° r°° , yl-\ \V2 

3 = 2^'J 1 k k dvidv2dv *-W 1 -i){yi~l)) Y 
= 0.0008144625656625044393912171285627219978 ■ • • (2) 

y = yi -2/3 

(yi + 2/2) (2/2 + 2/3) (yi + z/2 + 2/3) 

Focusing on x' 3 ': an d since this physical solution is known as a series expansion at 
w = (low or high temperature expansions), we will give all the connection matrices 
between this w = regular singular point and all the other regular singularities of 
the differential equation including the two new complex regular singularities |5] 
which are roots of 1 + 3u> + Aw 1 — 0. We will comment on the occurrence of the 
"ferromagnetic constant" 1% in the various blocks of the connection matrices. The 
decomposition of x^ i n the well-suited basis for each regular singular point allows 
us to find all the singular behavior of the physical solution. From these results, 
we will deduce the asymptotic behavior of the coefficients of the series expansion 
of x • These last problems are interesting, per se, for series expansions analysis 
of lattice statistical mechanics, since they correspond to subtle analysis of confluent 
singularities. Actually, we will see that even the last asymptotic evaluation problem is 
a (global) connection problem since the physical solution like x^ does not correspond 
to the obvious dominant singular behavior one might have imagined from the indicial 
equation. 

Focusing on the two new singularities, the roots of 1 + 3u> + Aw 2 = 0, we will 
show that the physical solution x * s not singular at these points. The factor of 
the logarithmic term, in the decomposition of x^ at these singular points, is known 
exactly and vanishes identically. 

Note that a fundamental concept to understand (the symmetries, the solutions of) 
these exact Fuchsian differential equations is the so-called differential Galois group |11| . 
Differential Galois groups have been calculated for simple enough second order, or even 
third order, ODE's (see for instance ^2). However, finding the differential Galois 
group of such higher order Fuchsian differential equations (order seven for x' 3 \ order 



Galois group for the Ising model 



4 



ten for x^) with eight regular singular points (for x^ 3 ') is no t an eas Y task J2| an d 
requires the computation of all the monodromy matrices associated with each (non 
apparent) regular singular point, considered in the same basis^. 

We will give the exact expression of all the monodromy matrices expressed in 
the same (w — 0) basis of solutions, these eight matrices being the generators of the 
differential Galois group, which will be given in a forthcoming publication |14| . 

This method can be generalized, mutatis mutandis, to the Fuchsian differential 
equation of x^- Here, we give the connection matrix between w — and, both, 
the ferromagnetic, and anti-ferromagnetic, critical points. The singular behavior is 
straightforwardly obtained with the asymptotic behavior of the series coefficients of 
the physical solution x • The monodromy matrices, expressed in the same basis of 
solutions are also obtained. 

The paper is organized as follows. We recall, in Section 2, some results on the 
Fuchsian differential equation satisfied by x j an d give a new factorization for the 
corresponding order seven differential operator, yielding the emergence of an order 
two, and an order three, differential operator (denoted Z 2 and Y 3 below). We 
give, in Section 3, the connection matrices matching the (series) solutions around the 
regular singular point w = and around all the other regular singular points. With 
these connection matrices we deduce the singularity behavior and the asymptotics 
on the physical solution of this ODE (Section 4). In Section 5, we deduce the 
exact expressions of the monodromy matrices expressed in the same basis. Section 6 
generalizes these results to the Fuchsian differential equation satisfied by x^ ■ Some 
physics implications of our results at scaling are discussed in Section 7. Our conclusion 
is given in Section 8. 



2. The order seven operator L-j 

Let us first recall, with the same notations as in [HUHji the seven linearly independent 
solutions given in [HIE] for the order seven differential operator L7, associated with§ 

x (3) . 

One finds two remarkable rational, and algebraic, solutions of the order seven 
differential equation associated with x^ 3 \ namely: 

S(L X ) = —0L-, 5(JV X ) = W * == = (3) 

1 - 4w (1 -4w) Vl - 16 w 2 

associated with the two order 1 differential operators given in [8]: 

Lx = — - 1 Nl = —- 2(1 + 2w) (4) 

dw w(l-4t»)' dw w (1 — 16 w 2 ) 

There is a solution behaving like w 3 , that we denote S3: 

S 3 = w 3 + 3 w 4 + 22 w 5 + 74 w e + 417 w 7 + 1465 w 8 

+ 7479 w 9 + 26839 w 10 +■■■ (5) 

and three solutions with logarithmic terms given by equation (17) in |Hj. Note the 
singled-out series expansion starting with w 9 , corresponding to the physical solution 

I These monodromy matrices are the generators of the monodromy group which identifies with the 
differential Galois group when there are no irregular singularities, and, thus, no Stokes matrices 1131 . 
§ x (n) is defined as x (n) = (1 - s 4 ) 1/4 /s ■ x (n) , for n odd. 
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S 9 = * 3) j W ^ = w 9 + 36™ 11 +4u> 12 + 884u> 13 + 196u> 14 + • • • (6) 
8 

The choice of this set of linearly independent solutions (and of these series) is, 
in fact, arbitrary since any linear combination of solutions is also a solution of the 
differential equation. Three of the above solutions are however singled out: the 
solutions S(L\) and S(N±) which are global (since they have closed expression), and 
the series Sg associated with the highest critical exponent in the indicial equation 
(w 9 +• • •), which has a unique (well-defined) expression and happens to correspond to 
the "physical" solution x^- Linear combinations, like S3 — a ■ Sg, are, at first sight, 
on the same footing. 

Nevertheless, introducing such a specific linear combination, B. Nickel|| has been 
able to show that the resulting series for the particular value a = 16 is, also, the 
solution of a linear differential equation of lower order, namely order four. With this 
result, the factorization scheme of L 7 becomes J : 

L 7 = Mi ■ Y 3 ■ Z 2 ■ Ni = B 3 ■ X 1 ■ Z 2 ■ Ni (7) 
= B3 • B 2 ■ 0\ ■ N\ = B3 ■ B 2 ■ T\ ■ L\ 

where the indices correspond to the order of the differential operators (B 3 , Y 3 are 
order three, B 2 , Z 2 order two, ...). The differential operators L7, Mi and Ti have 
been given in [S]. We give in Appendix A, the differential operators X\, Z 2 and I3. 
With these differential operators, all the factorizations Q can be found by left and 
right division. 

From these factorizations of L7, one can see that the general solution of the 
corresponding differential equation is the direct sum of the solution of Li and of the 
general solution of the differential operator Lg = Y3 ■ Z 2 ■ Ni. The operator L7 has 
the following decomposition: 

L 7 = L 6 ®Li. (8) 

We thus consider, from now on, the differential operator L 6 . 

The formal solutions of L e (at the singular point w = Q) show the occurrence 
of three Frobenius series and three solutions carrying logarithmic terms. With the 
factorizations Q, it is interesting to see which operator brings with it a singular 
behavior for a given regular singular point. Table 1 shows the critical exponents at 
each regular singular point for both differential operators Z 2 ■ Ni and Y3 ■ Z 2 ■ Ni. 
In the third and sixth column the number of independent solutions with logarithmic 
terms is shown. 

At the singular points w — 1, w = —1/2, and at the two roots wi, w 2 of 
1 + 3w + 4w 2 = 0, we remark that the solution carrying a logarithmic term is in 
fact a solution of Z 2 ■ Ni. Therefore, the three solutions of the differential operator 
Y3 ■ Z 2 ■ Ni, emerging from Y3, are analytical at the non physical singular points 
w = 1, w = —1/2, and at the quadratic roots of 1 + 3u> + 4u> 2 = 0. At the singular 
point w = 1/4, we also note that the differential operator Z 2 ■ Ni is responsible of 
the (1 — 4w) _1 behavior. We will then expect the "ferromagnetic constant" I3 to be 
localized in the blocks of the connection matrix corresponding to the solutions of the 
order three differential operator Z 2 ■ Ni at the point w — 1/4. 

|| We thank B. Nickel for kindly communicating this result. 

X The order four differential operator found by B. Nickel corresponds to B2 -Ti ■ Li = B2 ■ Ox ■ N\ = 
X 1 -Z 2 -N l . 
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Table 1: Critical exponents for each regular singular point for the differential 
operators Z 2 ■ N\ and Y% ■ Z 2 ■ N\. The columns TV show the number of solutions 
with logarithmic terms. The columns P show the maximum power of the logarithm 
occurring in the solutions. 

As far as explicit calculations are concerned, a well-suited basis necessary for 
cxplicitely writing connection matrices exists and can be described. Considering the 
order six operator Lq = Y3 • Z 2 ■ Ni , we construct the local solutions, sequentially, as 
the global solution of Ni then the two solutions coming from Z 2 - Ni, to which we add 
the three further solutions coming from Y 3 • Z 2 ■ N\ . We will use below this well-suited 
bases. 

3. Connection matrices for x' 3 ) 

Using a very simple method, let us show, in the case where one has an exact Fuchsian 
differential equation, that one can actually very simply, and very efficiently, obtain 
the connection matrices between two sets of series-solutions valid at two different 
points. The method consists in equating, at some matching points, the two sets of 
series corresponding, respectively, to expansions around w — and, for instance, 
w = 1/4. The matching point should be in the radius of convergence of both series. 
The singular points (i.e., w = and w = 1/4) should be neighbors, having no other 
singularity in between. Recall that the differential equation for x' 3 ^ has eight regular 
singular points, the point at infinity, five on the real axis and two (wi and w 2 ) on 
the upper and lower half plane each. At a given singular point w s , the solutions are 
obtained as series in the variable x, where x = w (resp. x = 1/w) for the point w s = 
(resp. w s — 00) and x = 1 — w/w s for the other regular singular points. We take the 
definition ln(x) = ln(— x) + iir for negative values of x which corresponds to matching 
points in the lower (resp. upper) half-plane for w > (resp. w < 0). 

The computation of the connection matrix should be more efficient when two 
"neighboring" singularities are, as far as possible, far away from the other singularities 
and, especially, when the test points chosen half-way are, as far as possible, far from 
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the other singularities, in order not to be "polluted" by the other singularities. We 
remark that one can calculate, in this way, just "neighboring" singularities: connection 
matrices of two singularities Wi , w r that are not "neighbors" should be deduced using 
some path of "neighboring" connection matrices: 

C(wi, w r ) = C(wi, w 2 ) ■ C{w 2l w 3 ) ■ ■ ■ C(w r _i, w r ) (9) 

This is the prescription we take for the singular points on the real axis and the 
singularity w\ lying in the upper half-plane. For the singularity w 2 lying in the lower 
half-plane, the connection matrix is calculated from: 

C(0, w 2 ) = C*(0, -1/4) • C*(-l/4, wi) = C*(0, wi) (10) 

where * denotes the complex conjugate. 

Let us remark that changing the variable w we are working with, to the 
more traditional s = sinh(2K) variable, or the usual high-temperature (resp. low 
temperature) variable t = tanh(if), or the variable r = (1/s — s)/2, modifies 
the distribution of singularities in the complex plane and their radii of convergence. 
However, the method can still be used. One can use that freedom in the choice of the 
expansion variable to actually improve the convergence of our calculations. 



3.1. Connecting solutions 

Let us first show, as an example, how we compute the connection matrix between two 
neighboring regular singular points (w = and w = 1/4) for order three differential 
operator Z 2 ■ N\. Around the singular point w = 0, the local solutions are two 
Frobenius series (one being the global solution S(N\)) and a series with a logarithmic 
term. The chosen basis is then (where x = w): 

S ( °\x) = 5(JVi)(aO, Si 0) (x) = [0,1,5,26,106,484,---], (11) 
Si°\x)= Si°\x)-In(x) + S^(x) (12) 

with: 

S30 ( x ) = [°, °> 0, 6, 26, 529/3, 2149/3, • • •] (13) 
where [ao, ai, »2, • • • ,] denotes the series ao + a\ x + a 2 x 2 + ■ ■ ■ There are three 
independent series s[°\ and 5^, since the operator Z 2 ■ Ni is of order three. 
Similarly, around w = 1/4, the local solutions read (with x = 1 — Aw and, where again, 

S^ 4 ^ is the global solution corresponding to operator Ni): 

S^ix) = SiNJix), (14) 

2 v ; x 4 96 64 55296 v ; 

S^\x)= si 1/4 \x)-Hx)+ S g /4 \x) (16) 

S 3o /4) ( x )= [3/8,-367/5760,-193/6720,-244483/6635520,- • •] (17) 

The series are defined around w — 0, and are convergent in a radius of 1/4, 
which corresponds to the nearest regular singular point (i.e., w = 1/4). Similarly, the 
solutions S' 1 ^ are convergent in the disk centered at w = 1/4 with same radius (i.e., 
1/4). Between the points w = and w = 1/4, there is a region where both sets of 



with: 
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solutions (S+ and S^ 1 ^') are convergent. This region corresponds to the common 
area between two disks centered respectively at w — 0, and w — 1/4, with the same 
radius 1/4. 

Connecting the local series-solutions at the regular singular points w = 0, and 
w = 1/4, amounts to finding the 3x3 matrix C(0, 1/4) such that 

S (0) = C(0,l/4)-S (1/4) (18) 



j(l/4) 



where (resp. S^ 1 / 4 -*) denotes the vector with entries (resp. S^ 4 ^). The 



.--(1/4) N 



solutions and Si 1 ' 11 ' are evaluated at three arbitrary points around a point x c 

belonging to both convergence disks of the series-solutions and S^ 4 ^. 

Equation (|18J) is thus a linear system of nine unknowns. The entries of the 
connection matrix C(0, 1/4) are obtained in floating point form with a large number 
of digits. These entries are "recognized" in symbolic form and matrix C(0, 1/4) then 
reads: 

1 



?(V4) 



(7(0,1/4) 



1 







9V3 

6477 



2 
3 

32 



ln(24)) 



64tt 





(19) 



The entries of this matrix are combinations of radicals, of powers of tt and 
logarithms of integers. Note that there is no straightforward manner to recognize 
numerical values such as the ones displayed above. However, it is possible, in a "tricky 
way" , to get rid of the logarithms of integers in the entries, and obtain as many zero 
entries as possible. This is shown, in the following, for this very example. 

The series, in the set of local solutions , are solutions of the differential 

equation (ODE) corresponding to the third order differential operator Z2 ■ N\ at 
the regular singular point w — 1/4. It is obvious that any linear combination of these 
series is also a solution of the differential equation. Consider the following combination 
instead of the third component in (|16(l : 



S 



(1/4), 



, 3 >{x) ^ (ln(x/2A) + 2/3)-Si 1/4 \x) +S&'*>(x) (20) 

By writing the argument of the logarithm as cc/24, there will be no logarithm in the 
connection matrix. Furthermore, by adding the second component of the basis to the 
third component with a factor of 2/3, the entry (2, 2) of the connection matrix will be 
canceled. The connection matrix then reads: 





1/4), 



(7(0,1/4) - 



1 



1 









37r%/3 
32 







(21) 



These tricks, based on well chosen linear combinations of the solutions, allow us 
to obtain as many zeroes as possible, and to get rid of the logarithms. They will be 
used in order to compute the connection matrix for Lq between the point w = and, 
respectively, w = 1/4, w = —1/4 and w = 00. 

The chosen well-suited basis of solutions, at each regular singular point calls for 
some comment. The factorization of the differential operator Lq being Y3 ■ Z2 ■ iVj, our 
method of producing the solutions, sequentially, allows one to determine from which 
differential operator a given solution emerges. Near the points w = 0, w = ±1/4, 
and w = 00, the third order differential operator Y3 brings three solutions (see Table 
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1), one Frobenius series, one solution with a log term, and one solution with a log 2 
term, denoted respectively §4, S5 and Sq . The solutions of the differential operator 
Y3 itself are of elliptic integral type (see Appendix B). These elliptic integrals behave 
around w — ±1/4 (resp. w — 00) like g(t) ■ ln(i/16) + f(t), with t = 1 — lQw 2 
(resp. t = l/16w 2 ), g(t) and f(t) being series with rational coefficients. One may 
then assume that the logarithmic term that appears in the solutions of Lq, inherited 
from Y3, will be of the form ln((l — 16w 2 )/16), near w — ±1/4, and of the form 
ln(l/256/w 2 ), near w = 00. The general form of combination for the fourth to sixth 
components of the well-suited basis will be: 

5 5 — ► S 5 + (ai - ln(c)) • S 4 (22) 

5 6 — > S 6 + 2 (a x - ln(c)) • S 5 + (ln(c) 2 - 1a x ln(c) + a 2 ) ■ S 4 

where c = 1, 8, 16 for the basis at, respectively, w = 0, w = ±1/4 and w = 00. The 
values of the parameters a\ and a 2 depend on each basis. 

Note that the argument in ln(x/24) in the series solutions of the differential 
operator Z 2 ■ N\ at w = 1/4 will be ln(x/4) and ln(x/24) at respectively w = 00 and 
w = 1. Similarly to I3, these arguments may come from the explicit solutions of Z 2 . 



3.2. Connection matrix between w = and w = 1/4 

The first three local solutions at w = are given by (|ll|l . 1)12(1. . and the fourth, 
fifth and sixth solutions read 

S ( °\x) = [0, 1,9,34,178,692,---], 

Si 0) (x) = S[°\x) ■ \n(x) + Si° \x) - sf(x)/4, 



Sr(x) = S™(x) - ln 2 (x) + 2 (S^(x) - S^(x)/i) ■ ln(x) 
? (o) 

7 50 



S${x)-S&\x)/2 + 25Si°\x)/16 



with: 



S$(x) = [0,0,0,-2,34,241/3, 



>50 



S^ ; (x) = [0, 0, 0, 0,-19/3, -7693/72, -575593/1800, • • •]. 

At the singular point w = 1/4, we make use of the combination (1221 which amounts 
to taking x/8 as argument of the logarithms in the fourth, fifth and sixth component. 
The parameters a\ and a 2 in 1221 are respectively 23/6 and 41/9. The first three local 
series at x = 1 — Aw are given in (UH), (El), jnj, (H, and the fourth, fifth and sixth 



read 



Si 1/4) (:r)= [1,-1/8, 3/16, 29/512, •••], (23) 
S^x) = (ln(x/8) + 23/6) • s| 1/4) (x) + sg /4 \x), 

S ( 6 1/4 \x) = (m 2 (x/8) + ^ ln(x/8) + 41/9^ • S| 1/4) (x) 



with: 



3 

+ 2 (ln(x/8) + 23/6) • S^ /4) (x) + sg /i} {x) 



S^ /4) (x) = [0,457/480,-2231/1680,-128969/184320,- • •] 
Seo /4) (x) = [0,-967/100,4312219/470400,595578701/116121600,- 
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Connecting both solutions amounts to solving a linear system of 36 unknowns (the 
entries of the connection matrix) . We have been able to recognize these entries which 
are obtained in floating point form with a large number of digits. The connection 
matrix C(0, 1/4) for the order six differential operator Lq reads: 

£7(0, 1/4) = (24) 



±11 

16 



1 














o ■ 


1 





9^ 
64tt 














3?r\/3 
32 














5 


i - 2 • 7+ 

3 Z J 3 


3V3 
64tt 








1 


5 
4 


371"^ 

32 


45^ 
256tt 





1 

32 





o 2 

Ztt 
3 


64 


225\/3 3tt\/3 
1024tt 64 


n 1 
64 









Some comments on how these entries have been "recognized" will be given below. 
Let us remark that, once the entries of the connection matrix have been obtained, a 
further change of basis can be made to get it as "simple" as possible. 



3. 3. Connection matrices between w = and the other regular singular points 

The chosen basis of solutions and the connection matrices between w — (high or low 
temperature) and, respectively, the anti-ferromagnetic point w = — 1/4 and the point 
w = oo (corresponding to s = ±i) are given in Appendix C. 

The chosen basis, used for the regular singular points w — 1,-1/2 and 
l + 3u> + 4w 2 = 0, are given in Appendix D together with the corresponding connection 
matrices with the point w = 0. Many entries are "recognized" and, in particular, those 
required to find the singular behavior of the physical solution. They correspond to the 
third column of matrices given in Appendix D. 

The connection matrix between each pair of neighboring singular points is 
computed with the well defined procedure described above. The connection matrix 
between w = 0, and a non neighbor singular point, is computed using For instance, 
£7(0, 1) is computed from £7(0, 1/4) and £7(1/4, 1) as £7(0, 1) = £7(0, 1/4) • £7(1/4, 1) 
which says that the solutions defined at w = 1/4 connected to the solutions defined 
at w = 0, are also the solutions that are connected to the solutions defined at w = 1. 

To be more confident of this prescription, let us underline that the connection 
matrices C(0, 1) and C(0, — 1/2), deduced from ©, will be used below to confirm 
known dominant singular behavior of an d find the subdominant behavior. 



3.4- Comments and remarks 

The connection matrices between w = and the other singular points are structured 
in blocks. The latter, due to the factorization of the differential operators and to 
the sequential building of the solutions, are easily recognized. The block (1,2,3) x 
(1,2,3) is associated with the third order differential operator ■ N\. The block 
(4, 5, 6) x (4, 5, 6) represents the connection between the solutions (at both w = and 
the other singular points being considered) of Lq that are not solutions of Z2 ■ N±. 
The "ferromagnetic constant" 1^ appears in the connection matrix between w = 
and w = 1/4, as mentioned earlier, in the block (1,2,3) x (1,2,3) at the column 
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corresponding to the S 2 (see l|15[l) solution of the third order differential operator 
Z 2 -N 1 . 

To compute the connection matrix, we have used the differential operator Lg 
which has a unique factorization. If, instead, we consider the differential operator L7, 
the next solution (around w = 0), that comes from Mi, will be the series (jHJ and 
will correspond to x^- This seventh solution is expressed as a linear combination of 
the already existing components and of the solution of the differential operator L\. 
We can then choose to add the latter as the seventh solution. The connection matrix 
will have a 1 at the entry (7, 7) and zero elsewhere on the seventh line (and column) , 
since the solution of the differential operator L\ is global. By considering another 
factorization of L7, we will get the same structure with an obvious relabelling. 

Let us make a few computational remarks on the calculation of these connection 
matrices. At the matching of the series-solutions for which 1500 coefficients§ are 
generated from homogeneous and non-homogeneous recurrences, the entries of the 
matrix are computed with 800 digits for all the singular points. The numbers that 
come in floating form are "recognized" as powers of 7r, radicals and rational numbers, 
and are in agreement up to 400 digits | for the connection between the solutions at 
w = and w — ±1/4, and up to 100 digits for the connection involving other singular 
points like, w = 1. This fact is related to the convergence rate of the series at the 
(midway) chosen matching points. For instance, between w = and w = 1/4, the 
matching points near w = 1/8 are such that both series (at w = 0, and w = 1/4), 
which have the same radius of convergence, will be faithfully reproduced with the 
number of terms used in the series. The matching of the solutions between w = 1/4, 
and w = 1, will then require more terms to fulfill the same accuracy than in the 
(w = 0)-(w = 1/4) situation. This is due to the fact that, at w = 1, the convergence 
radius being 3/4, the matching points, which should be in the common region of both 
disks, are closer to w — 1/4 than to w = 1. As a general rule, the matching points are 
chosen around the middle of the segment in the common area between the convergence 
disks of the two regular singular points for which the connection matrix is computed. 

The difficulty in finding "non-local" connection matrices is rooted in the 
recognition of the entries. We have given the connection matrix between w = and 
w = 1/4 with entries fully recognized (apart from 1%) to show that the method actually 
works and is efficient. For the matrices concerning the connection between w = and 
the other singular points, we have concentrated our effort on the entries that will show 
up in the physical solution. We should note that there is no reason to expect the other 
(not yet recognized) entries to be "simply" combinations of it's, log's and radicals. 
These entries are probably valuations of holonomic functions. This was clearly seen 
in numerous examples we tackled of various differential equations (of order two and 
three) with known solutions of hypergeometric type. The recognition process used the 
fact that we actually found the explicit solutions of differential operator Y3 and, thus, 
knew how the numerical logarithms can be tackled. These were "absorbed" in the 
basis. We know, on the other hand, that the problem is roomed with hypergeometric 
functions. We then expect some it's to be present. For the entries consisting of simple 
product expression, recognizing the number amounts to performing simple arithmetic 
operations. Note that considering the inverse of the connection matrix, some entries 
also show up as simple rationals. The combination where 7r's, radicals and rationals 

§ For some checks, 3000 terms have been generated. 

|| Let us note that the "ferromagnetic constant" has been obtained up to more than 400 exact 
digits. 
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appear additively comes from looking to, for instance, at the determinant of the 
matrices, or block matrices, which happen to be easily recognizable (in fact rational 
or quadratic numbers for the roots of 1 + 3w + Aw 2 = 0). 

Another remark is the following. We first obtained the connection matrix l|24|l in 
some general basis. The matrix had more non zero entries compared to (|24|l involving 
powers of 7r, radicals and also ln(3) and powers of ln(2). The well-suited basis we chose 
has "evacuated" all these log's in the entries of the matrix, lessening the recognition- 
process effort. But, of course, all these logs will reappear in the final result such as 
the singular behavior of the physical solution as next sections will show. 



4. The physical solution x' 3 ^ and its singular behavior 

The calculations of connection matrices are obtained straightforwardly from the well- 
defined numerical process described in Section Having N singularities, one needs 
N — 1 such connection matrices in order to find the correspondence between all these 
well-suited bases of series-solutions. 

Let us focus on some particular entries of these various connection matrices, 
namely the entries corresponding to the decomposition of x^ m terms of the various 
well-suited bases associated with each singularity. We have used the fact that the 
physical solution (corresponding to X ) decomposes as the solution of differential 
operator L\, S(Li) (which is X /%) and the physical solution of the operator Lq 
denoted $ 6 (w) 00: 

X (3) H - ^X (1) + *e(«0 

Furthermore, our well-suited basis of solutions at the singular point w — 0, does not 
contain, as a component, the physical solution ^e(w) which is given in terms of the 
previously considered components as: 

*.(»)= jSF-^sW-lsF (25) 

This physical solution can now be easily obtained from the connection matrices 
between w — and any regular singular point, that we denote w = w s (with x = w, 
x = 1/w for respectively w = and w = oo and x = 1 — w / w s , otherwise) as: 

6 - 

= E (3 - u c(0 ' Wshl ~ I c{0 > Ws)4j ) ' s i Ws) 

For instance, at the ferromagnetic critical point, this physical solution ^e(x) can easily 
be deduced from J23J, and written as: 

&„( x \ - _i (I -21+)- S (1/4) — 9 (1/4) 



4 1/4) and 4 1/4) are known from their series expansion HSJ, This equation, 

giving the full expansion of x at w = 1/4, can hardly be obtained directly from 
the integrals defining x < " 3 H w )- O ne nas similar expansions for all the other singular 
points. 
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4-1. Singular behavior of 

Knowing the behavior of solutions S^"^ near each regular singular point, it is 
straightforward to get the singular behavior at those points for the physical solution 
<fr 6 (and thus x^)- 

Considering the critical behavior of x*- 3 ' near the ferromagnetic critical point 
w = 1/4, and denoting x = 1 — Aw, the singular part of the "physical" solution x^ 3 ' 
reads: 

X {3) (singular, 1/4) = ± £ - ^ s| 1/4) • ln 2 ^) (26) 

+ ^((3 1n(2)-|).5| 1 / 4 )-^ /4) )-ln(x) 

where l£ is actually the "ferromagnetic constant" J3J), and S^^ the series defined 
in the well-suited basis l)23|) at w — 1/4. The results agree with previous results of 
B. Nickel, but the correction terms are newt, m particular the term 31n(2)/32/7r 2 in 
H2(jf> . In terms of the t = (1/s — s)/2 variable introduced in |H>1 1151 [TB] . the singular 
part H26fl reads: 

-fin,- i It ln 2 (r) /. 23 \ ln(r) 

X (3) (singular, r ^ 0) ~ ± - ^ + (ln(2) - - J • -±± + • • • 

Near the antiferromagnetic critical point w — —1/4, x' 3 ) behaves as: 

X< 3 > (singular, -1/4) = - J-j ^ 1/4) ■ ln 2 (x) (27) 

-^(3(2-ln(2))-Si- 1/4) +4o 1/4) )'M-) 

At the non-physical singularities w = 1 and w = —1/2 the physical solution 
behaves, respectively, like: 

X™ (singular, 1) = ■ ■ ln(x) (28) 

and 

X (3) (singular, -1/2) = -f^ • S^ 1/2) • ln(x) (29) 

confirming Nickel's calculations given in |7|. 

At the point w = oo, corresponding to the non physical singularities s — ±i, the 
singular behavior reads: 

X (3) (singular, oo) = - -L S M ■ In 2 (*) (30) 
lbn z 

- ± ((4 - 27rl ) • - (5 + 41n(2) + i |) • S<°°> + • ln(z) 

At the new singularities found in |SJ, namely the roots of 1 + 3w + 4w 2 = 0, 
which are regular singular points of the differential equation, the singular part of the 
physical solution reads, at first sight: 

X (3) (singular, wi) = - — (a 23 + 3a 43 ) • S ( 2 Wl) ■ ln(x) 

The entries 023 and 043 (see the connection matrix for these points in Appendix D) 
are however such that 023 + 3a 4 3 = 0. The physical solution is thus, not singular, at 
the newly found quadratic singularities, confirming our conclusion given in 9 from 
series analysis. 

I These results have also been found by B. Nickel (private communication). 
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4-2. Asymptotic series analysis 

As the physical solution is given as a series around w = 0, the coefficients of 
the latter are controlled by the nearest singular points (i.e. w = ±1/4). Since the 
singular parts at the ferromagnetic and anti- ferromagnetic critical points l|26|l . (|27ll 
are obtained, it is straightforward to deduce the behavior of the coefficients of series 
10 for large values of n. Standard study of the asymptotic behavior of the coefficients 
via their linear recursion relation can be used (see |17|). For our purpose, we use the 
following identity for ln 2 (l — x) (where x stands for x = Aw): 



ln 2 (l - a;) = b(n) ■ x n , where: 

n— 2 

n— 1 , n 

E-7^7T= -(nn) + l) (3D 



^-^ i (n — i) 

i=i y ' 



where 7 = 0.57721566 • • • denotes Euler's constant, and 'J denotes the logarithmic 
derivative of the T function. Recalling the asymptotic expansion of ^f(n) up to l/n 2 
for large values of n, one obtains: 

With the same manipulations of ln 2 (l + x), and inserting in 1261) . (|27|l . one obtains 
the asymptotic form of coefficients of /8w 9 as: 

0-15 c(n) _ I 3 + 1/1 , , /Mn) , h 1 



1 rf5 + .(-l)-U + 



167r 2 V2 V ' J \ n n 2n 2 



16tt 2 V12 
where 61 = 7 + 31n(2). 

It is this parity effect in the asymptotic behavior of the coefficients that we 
saw, numerically, (see equations (33) in [5]) where we obtained, around n ~ 500, 
c(n) ~ 13.5 x 4™ for n even and c(n) ~ 11 x 4™ for n odd. For very large values of 
n, the asymptotic value of the coefficient c{n)/A n is thus 2 14 ■ /g 1- ~ 13.34415467 • • •. 

5. Monodromy matrices for 

5.1. Sketching the differential Galois group of 

As a consequence of the direct sum (jHJ , the differential Galois group of L7 reduces (up 
to a product by C) to the differential Galois group of Lq. From the factorization of Lq, 
one can immediately deduce that the differential Galois group of L§ is the semi-direct 
product of the differential Galois group of I3, of the differential Galois group of Z2 
and of the differential Galois group of N\ (namely C). 

In some "well-suited global basis" of solutions, the form of the 6x6 matrices 
representing the differential Galois group of L 6 , reads: 



A 
H B 



with 



b 
h g 



where the 2x2 matrix g, and 3x3 matrix B correspond, respectively, to the 
differential Galois group of Z2 and Y3. The 3x3 matrix A is associated with the 
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differential Galois group of Z 2 ■ iVj., and the 3x3 matrix H corresponds to the fact 
that we have a semi-direct product of the differential Galois group of Y 3 and Z 2 ■ N% 
in L 6 = L 3 ■ Z 2 N v 

Many papers (for instance |181 1191 121)1 1211 122j ) describe how to calculate the 
differential Galois groups of order 2 and order 3 differential operators. The differential 
Galois group of Lj will be deduced in a forthcoming publication |14| . 

To go beyond this sketchy description of the differential Galois group, one needs 
to calculate specific elements like the monodromy matrices expressed in a common 
basis. 



5.2. Monodromy matrices rewritten in the w = basis 



Having the connection matrices between w = and each singularity, the local 
monodromy matrices expressed in their own well-suited basis of (series) solutions, 
can be rewritten in a unique global basis valid for all singularities. This will allow us, 
in a second step, to calculate their products and thus generate the differential Galois 
group. Let us define the 2x2 and 3x3 matrices 



A 



1 



B 



1 

Q 

n 2 



o 
i 

2Q 



(32) 



where denotes 2 i tt and corresponds to the translation of the logarithm when 
performing a complete rotation around the regular singular point: ln(w) — > ln(w) + J7. 

The expression of the local monodromy matrix around each regular singular point 
w s in its own well-suited basis of (series) solutions reads: 

e 

l(w s ) = C (33) 
D 

where, e and the 2x2 blocks C, and 3x3 blocks D, are such that: 

w = 0, w = oo -> e = +1, C = A, D = B 

w = 1/4, -> e = — 1, C = A, D = B 

to = -1/4, -> e = -1, C = Id, D = B 

w = 1, -1/2, -3/8 ± i Vf/8, -> e = +1, C = A, D = Id 
The monodromy matrix around any singularity w — w s expressed in terms of the 



(w = 0) well-suited basis, and denoted M w= o(w„), reads: 

M„=oK) = C{0,w s ) ■ l(w s )(Q) ■ C-i&Ws). 



(34) 



In order to keep track of the n corresponding to the translation of the logarithm in 
the local monodromy matrix /(iy s )(f2), and the 7r's occurring in the expression of the 
entries of the (quite involved) connection matrix C(0, w s ), we will denote the latter 
by a = 2 in. 

Let us focus on the singular point w = 1. Its monodromy matrix, expressed in 
terms of the w — well-suited basis, is given by l|34l) with w s = 1, and where the 
connection matrix C(0, 1), matching the (w = 1) well-suited basis together with the 
(w = 0) well-suited basis, is a "quite involved" matrix given in Appendix D, with 
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entries depending on 7r's and on a set of 15 constants, not yet recognized in closed 
form. The monodromy M w= q(1) can finally be written as a function of only a and 
ft: 



8 a 2 ■ M w=0 (l)(a,Q.) = 
8 a 2 

-48 aft 



-1008 a ft 
12 a (5 + 16 a) ft 



(35) 





8 a 2 








—48 ^ 
8 a 2 
-1008 ft 
12 (5 + 16 a) ft 
- (75 + 44 a 2 ) ft 






























a 2 











8 a 2 











8 a 2 



-a (75 + 44 a 2 ) ft 

Let us give one more example corresponding to the new quadratic singularities 
1 +3w +4w 2 = 0. The monodromy matrix around one of the quadratic singularities 
w = W\, expressed in terms of the (w — 0) well-suited basis, after the conjugation 
(|21Jl, reads: 



8a 2 -M w=Q (w 1 )(a, ft) 



A 
B C 



(36) 



with: 



.4 
B 



8 a 2 
48 aft 
16 ft a 2 
-16aft 











8 a (a + 6 ft) -144 ft 

16 ft a 2 -8 a (-a + 6 ft) 

-16 a ft 48 ft 

4a (4a- 15) ft 4a (4a- 15) ft -12 (4a- 15) ft 

aa aa —3a 

with a = (-40 a + 12a 2 + 75) ft and [ C ] = 8a 2 • Id(3 x 3). 

One can actually verify that the monodromy matrix around the other quadratic 
singularity w — u>2 (w2 is complex conjugate of u>i), expressed in terms of the 
(w = 0)-well suited basis, actually identifies with (|36|l where a has been changed into 
—a. 

We have totally similar results for all the other (regular) singularities. The 
expression of the other monodromy matrices M w= o(w s ), around the other (regular) 
singular points w — w s , are displayed in Appendix E. 

We saw that the connection matrices depend on and on "still not yet 
recognized" (probably transcendantal) numbers, like X42 and yn (for the connection 
matrix between w — and w = 00). Rewriting a monodromy matrix in a 
unique (global) basis like the w = basis, amounts to performing conjugation, like 
H34J) . of simple (local) monodromy matrices depending only on ft, by these quite 
involved connection matrices. As a consequence, one does expect, at first sight, these 
monodromy matrices, rewritten in the unique w = basis, to be dependent on the still 
unknown numbers. For instance one certainly expects the monodromy matrix around 
w = 1/4 (see Appendix E) to be expressed in terms of the transcendental number 1^ , 
or the monodromy matrix l|35|) to depend on 15 parameters. It is worth noting that 
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all these matrices M(w s ), expressed in the same (w = 0) well-suited basis, turn out 
to be quite simple matrices where the entries are actually rational expressions, with 
integer coefficients of a and SI. Section (|5.3[1 gives some hints on why this is so. 

The introduction of the two parameters a and SI is a nice "trick" to track the 7r's 
coming from the connection matrices versus the 7r's coming from the local monodromy 
matrices. However, one should keep in mind that a is not independent of SI: the 
"true" monodromy matrices are such that a = SI (fi being equal to 2 in). Let us 
denote these "true" monodromy matrices by Mi, i 



Mi 
M 3 
M 5 



M u 
M„ 



M 7 = M u 



Ma 



The matrices M 2 
share the same Jordan block form, 
respectively: 

" -1 



= 1) • * " ) 8: 

(oo)(n, n), m 2 = m w=0 (i){si, si), (37) 
(i/4)(o, si), m 4 = M w=0 ( Wl )(n, si), 

o(-l/2)(n, SI), M 6 = M w=0 (-l/4)(Sl, SI), 
o (0)(fi, fi), M 8 = M w=Q {w 2 )(Sl, SI) 

M 5 , M 8 , and respectively the matrices Mi and M-j, 
The Jordan block forms for M3 and Mg read 







110 

110 

1 

1 1 

1 



-1 

110 

110 

1 

1 

1 



These matrices Mi are the generators of a 6x6 matrix representation of the 
differential Galois group of the Fuchsian differential equation corresponding to Lq. 
Any element of the differential Galois group is of the form: 

(38) 



M ni • M" 2 • M™ 3 

JW P(1) 1V1 P(2) IVI P(3) 



M P{4) 1V1 P{b) 1U P(6) M P(T) 1U P(S) 



where P denotes an arbitrary permutation of eight elements and where the n^'s are 
positive or negative integers. This looks, at first sight, like an infinite discrete group, 
but the closure of this infinite set of matrices can be quite large continuous groups 
like semi-direct products of SL(2,C) with SL(3,C), ... 

Our "global" (800 digits, 1500 terms) calculations yield quite involved exact 
connection matrices. With such large and involved computer calculations there is 
always a risk of a subtle mistake or misprint. At this stage, and in order to be "even 
more confident" in our results, let us recall that the monodromy matrices must satisfy 
one matrix relation which will be an extremely severe non-trivial check on the validity 
of these eight matrices Mi, or more precisely their (a, Si) extensions. Actually it is 
known (see for instance Proposition 2.1.5 in |23p. that the monodromy group§ of a 
linear differential equation (with r regular singular points) is generated by a set of 
matrices 71, 72, ■ • • , j r that satisfy 71 ■ 72 • ■ • 7r = Id, where Id denotes the identity 
matrix. The constraint that "some" product of all these matrices should be equal to 
the identity matrix, looks quite simple, but is, in fact, "undermined" by subtleties of 
complex analysis on how connection matrices between non neighboring singular points 
should be computed. The fact that the prescription H9I10|I has given no contradictory 

§ Which identifies in our Fuchsian case to the differential Galois group. 
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results on the singular behavior may be an argument that our Mj's are not "too 
far" from these "elementary" 7j's. In other words, one of the products l|38|) must be 
equal to the identity matrix for some set of n^'s and for some permutation P. With 
the particular choice lpT7)l of ordering of the eight singularities, this product, actually 
reads: 

Mi ■ M 2 ■ M 3 ■ M 4 • M 5 ■ M e ■ M 7 ■ M 8 = Id (39) 

Of course, from this relation, one also has seven other relations deduced by cyclic 
permutations. It is important to note that these relations (|39|) are not verified by 
extensions like <|35|) . i|36|) depending on two independant parameters a and f2, of the 
monodromy matrices Mj. If one imposes relations (|39|l for the (a, f2) extensions of 
the Mj's, one will find that, necessarily, a has to be equal to O, but (of course||) one 
will find that these matrix identities are verified for any value of 57, not necessarily 
equal to 2 i it. 

5.3. Comments 

The entries of the connection matrices have been seen to be expressed as various 
polynomials, or algebraic combinations of power of 7r, ln(2), ln(iV) (jV integer), 
algebraic numbers, etc., and more "involved" transcendental numbers like J2J). On 
the other hand, the monodromy matrices M w= q(w s ), expressed in the same (w = 0) 
well-suited bases, have entries which are rational expressions with integer coefficients 
of a and f2. To get some hint as to how this occurs, let us consider, for instance, the 
regular singular point w = 1 . The local monodromy matrix is almost the unity matrix 
(only one solution with log) with elements: 

= % + n • *3 S j2 (40) 

The product H34J) giving the global monodromy matrix will be given by 

M w=0 (% = +ft-C(0,l) l3 -C- 1 (0,l) 2 , (41) 

where one can see that only the third column of C(0, 1) and the second row of its 
inverse will contribute. These entries have been "recognized" (see Appendix D). 

Let us assume that there is another solution with a log term (this is not so, see 
Table 1). An entry (for instance /(l)6s) of the local monodromy matrix changes from 
zero to f2. In this case equation Q41JI becomes: 

Mv^Ma = Sij + n-c(o,i) t3 c- 1 (o,i) 23 + ft-c(o,i), 6 c- 1 (o,i) 5 , 

The entries C(0, 1)^6 and C _1 (0, 1)5^ will appear in the global monodromy matrix. In 
fact, changing the entry Z(l)65 from zero to means that a formal solution will exhibit 
log's, and this will correspond to the entries C(0, As a pratical rule, we found 
that such entries (corresponding to solutions with log's) can be easily "recognized" in 
contrast with the entries corresponding to Frobenius series which will be canceled by 
the zero entries of 1(1). The entries corresponding to Frobenius series are probably 
valuations of holonomic functions. 

Let us now assume (for the actual situation) that the whole column C(0, l)j3 has 
unknown entries. Recalling the fact that the product of the monodromy matrices, 
expressed in the same basis, should be equal to the identity matrix |23| (this is what 

I A matrix identity like 1391 yields a set of polynomial (with integer coefficients) relations on 
Q = 2 in. The number n being transcendental it is not a solution of a polynomial with integer 
coefficients. These polynomial relations have, thus, to be polynomial identities valid for any Q. 
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we found for our eight matrices M i} see H39[0 . one then expects the "not yet guessed 
constants" (i.e., the column C(0, 1)^) to be given by a non linear system of equations. 
This is indeed what occurs for this example, and we recover that way the entries given 
for this case in Appendix D. 

A last remark is the following. Right now, we have considered all the matrices 
(connection and therefore monodromy matrices expressed in a unique basis) with 
respect to the (w — 0) well-suited basis of solutions. This is motivated by the 
physical solution \^ which is known as series around w = 0. In fact, we can switch 
to another w = w well-suited basis of solutions. This amounts to considering the 
connection C(w, w s ) — C" 1 (0, w) ■ C(0, w s ). For instance, we have actually performed 
the same calculations for the (w = 1/4) basis of series solutions. We have calculated 
all the connection matrices from the (w — 1/4) basis to the other singular point basis 
series solutions, and deduced the exact expressions of the corresponding monodromy 
matrices now expressed in the same (w — 1/4) basis of series solutions. It is worth 
noting that we get, this time, for the monodromy M w= i/^{w s ) around singular point 
w s and expressed in the (w = 1/4) basis, a matrix whose entries depend rationally 
on a, f2, but, this time, also (except for the monodromy matrix at w — 1) on the 
"ferromagnetic constant" Jg . One verifies that the product of these monodromy 
matrices in the same order as (|39|) . is actually equal to the identity matrix when 
a = Q, the matrix identity being valid for any value of a = f2 (equal or not to 2 in), 
and for any value of 7g (equal or not to its actual value given in (0)). 

We have similar results for the monodromy matrices around singular point w s , 
expressed in the (w = oo) basis, but, now, the monodromy matrices M w — 00 {w s ) 
depend on a, fl, and, this time, on the (not yet recognized) constants yn and £42. 
Again, the product of these monodromy matrices in the same order as <|39|) . is actually 
equal to the identity matrix when a = Q, the matricial identity being valid for any 
value of a — (equal or not to 2 in) and for any values of j/41 and X42 (equal, or 
not, to their actual values given in Appendix C). 

6. Mutatis mutandis: Connection matrices and singular behavior for 

6.1. Connection matrices 

The Fuchsian differential equation forj x^ 4 \ the four-particle contribution to the 
susceptibility, is given in jl()| . The order ten differential operator £10 associated with 
this differential equation has 36 (equivalent up to isomorphisms) factorizations (see 
Appendix F in JO]). Consider, for instance, two of these factorizations: 



The notations are the same as those in JU], the M operators are of order four, the 
N and L operators are respectively of order two and one. The two factorizations 
above mean that £10 is a direct sum of an order eight differential operator, £g 
= AI2 ■ L25 ■ L12 ■ L3 ■ Lq and of the order two differential operator jVq (which, see |10j . 



£10 



Ns ■ M2 ■ L25 ■ L12 ■ L3 ■ Lq 

Mi ■ L24 ■ Li 3 ■ Ln ■ Lu ■ N 



(42) 



has remarkably as solution): 



£10 = £ 8 © iV 

I x' n ) is defined as x'™' = (1 — s~ 4 ) 1//4 • x'"'; f° r n even. 



(43) 
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As was the case for , it is thus sufficient to consider the differential operator £ 8 for 
which a general form of 8 x 8 matrices, representing Qal(£ s ), the differential Galois 
group of £ 8 , is deduced: 

L 

G M 

G, M and L are 4x4 matrices, the latter being lower triangular. Recall that Cs has 
four known global solutions (see ^H] and below). 

Similarly to the calculation on j£ , we can J f°r instance, calculate connection 
matrices associated with the correspondence between the series near x — 16 w 2 =0 
(high temperature) with the series near x — 16 w 2 =1 (ferromagnetic and 
antifcrromagnetic critical point), and find how the "physical solution" can be 
decomposed on the various well-suited bases around each singular point (physical or 
non-physical) of the order ten Fuchsian differential equation. 

We use the factorization i|4*2|) to construct the basis of solutions, sequentially, as 
the four solutions corresponding to the differential operator L 25 ■ £12 ■ L 3 ■ L that 
we call respectively Si, S2, S3 and S4. To these solutions, we add the four solutions 
coming from C% and inherited from the differential operator M 2 , that we call S5, • • •, 
5*8. Here, again, an optimal choice of basis is made in order to have as many zeroes 
as possible in the connection matrix with as "simple" entries as possible. The basis 
of solutions at x = and x = 1 (with respectively t = x and t = 1 — x) have similar 
forms and read: 

Si(t) =1, S 2 (t)= eq.(33) in [TU], 

S 3 {t) = eq.(32) in JTHj, S 4 (t) = eq.(43) in [Ej, 

S 5 (t) = see below, S 6 {t) = S 5 (ln(t/16) + a x ) + S 6a 

S 7 (t) = (\n(t/16) 2 + 2ai ln(t/16) + a 2 ) ■ S 5 

+ 2S 60 (ln(t/16) + 01) + 5to 
S s (t) = (ln(Vl6) 3 + 3 ai ln(t/16) 2 + 3a 2 ln(t/16) + 03) • S 5 

+ 3 (ln(t/16) 2 +2a 1 ln(i/16) + o 2 ) ■ S 6 o 
+ 3 (ln(t/16) +ai) ■ S 70 + S 80 

where the constants ai, a 2 and 03 and the series read, near x = 

ai = 79/60, a 2 = -751/1800, a 3 = -10619/375, 

Si°\t) = [0,0,1,45/32,425/256,945/512,- • •], 
S$(t) = [0,2/3,0,2353/13440,121619/322560, •••], 
Sm(t) = [8,-119/45,0,-560333/1411200,- • •], 
S$(t) = [0,0,0,0,-127639044817/85349376000, •••] 

and, near x = 1 : 

ai = 35/6, a 2 = 107/9, a 3 = -1051745657/749700 
S^\t) = [1,-1/4,-7/64,-45/256,-3385/16384, •••], 
S$(t) = [0,7/120,-3809/13440,42401/16120,9271027/18923520,- • ■], 
Srnit) = [°. 1099/75, 741847/78400, 218499331/101606400, • • •], 
S$ (t) = [0, 0, 0, -37462660457/592220160, • • •] 
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The connection matrix between x = and x 

1 

C(M) A B 



1 comes out as: 



(44) 



where 1 denotes the 4x4 identity matrix and denotes the 4x4 zero matrix. The 



4x4 identity matrix corresponds to the fact the four solutions Si, 
solutions. The two lower 4x4 blocks read: 



S4 are global 



A = 



051 


052 


5 
2 


a 5 4 





h 





1 

32 


a n 










asi 


-7T 3 


083 


a 84 



with 



a 71 = -T 



2422 



225 



5rr 2 
~6~ 



B 



2422 







_i_ 

2n 



1 

27T 3 







225 



a 84 = — 777: 



7T 

32 



1211 



600 



073 = 

The "not yet recognized" entries of this matrix read: 

<z 5 i ~ - 17.882936774520, a 52 ~ 7.767669067696, 
a 54 ~ 0.530951641617, a 81 ~ -92.773462923758, 
a 83 ~ 77.887072991056 

Here again, the block structure of the connection matrix relies on the factorization 
of Cg and on the "sequential" building of the solutions. The block matrix B represents, 
specifically, the connection between the solutions inherited from M% at both points 
x = and x = 1. This fourth order differential operator M% in £ 8 (corresponding to 
X^) is structurally very similar (see the remark at end of Appendix B) to operator Y3 
in Lq (x^ ) ■ Similarly to x^ case, a ferromagnetic (and anti-ferromagnetic) constant 
(see (|48t below) is localized at the fifth line. 

We have also computed the connection matrices§ (not given here) between the 
solutions at x — and respectively x — 4 (corresponding to Nickel's non-physical 
singularities) and x = 00 (corresponding to the non-physical singularities s = ±i). 
Denoting by M x= o(0), M x= o(l), M x= q(4) and M x= q(oo), the monodromy matrices 
expressed in the same x = well-suited basis obtained with similar conjugation like 
l|34(l . one obtains: 

M a=0 (oo) • M x=0 (4) ■ M x= o(l) • M X=0 (Q) = Id (45) 

This identity is valid irrespective of the still unknown constants. 



6.2. Singular behavior of x^ 

The particular physical solution corresponding to x*- 4 '' = X^ /3 + $s (see ^01) is given, 
in terms of the basis chosen at the point x = 0, by: 

$ 8 = ± -(5Sf -5Sf -2Sf) (46) 

At the ferromagnetic, and anti- ferromagnetic, critical point x = 1, the solution can be 
deduced from the above connection matrix and reads: 

$8 " ~384 ( 2a 5i-5)-S 1 -^"^ ~ W2 -S, +^3-^ 8 

§ The matching points are taken in the lower half-plane of the variable x. 
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Here again, the above decomposition corresponds to an expansion at the point x = 1 
of the triple integral defining x^ . 

From this solution, the singular part of x^ reads (with t — 1 — x): 

(singular, 1) = ^ + S P ■ In 3 (t) (47) 



-(im(2)-!).s# + ^s£>).m(t) 



+ ^2^(1/2, -1/2; 2; i)-ln(t) 
The constant £Q 77/ reads, in terms of the "not yet recognized" numbers 052, 054: 

II = — + ^| - ~ 0.0000254485110658 • • • (48) 

O07T Izo 24 

The first term at the right-hand-side of (|48J) comes from \^ j as well as the last term 
in (gZj). 

Similarly, the singular behavior of the physical solution at the other singular 
points can easily be obtained from the corresponding connection matrices (not given 
here). At the singular point x = 4, the physical solution behaves like (with t = 4 — x): 

X (4) (singular, 4) = - ' (l + -t + —t 2 + ■ ■ ■) (49) 

A v B ' ; 2 10 • 3 2 • 5005 V 4 272 J y ' 

confirming the calculations in 0. 

The singular behavior of x^ at the singular point x = 00 reads (with t = 1/x): 

X (4) (singular, 00) = -20 i ■ r 1/2 ■ (a + 3A X ■ \n(t) (50) 
+ 3(( 0l - 41n(2)) • S? + Sgg) • ln 2 (t) + Sf • ln 3 (i)) 



n=0 



with 



At = |(2Jf - 1) ■ S£ + (l61n 2 (2) - 801 ln(2) + a 2 ) ■ S 5 ° 



_>oc 
'70 



+ 2(oi-41n(2))- S™ + 3Sj 
A = 2n 3 (i 52 + i ^(2K - 1)) S 2 °°/5 

- f-i|(2X-l)+2(5 + 2r 53 ) N )7r 3 53 00 /5 



(64 In 3 (2) - 48ai In 2 (2) + 12a 2 ln(2) - a 3 ) • S5 
(2if - 1) • 5|§ + 3 (l6 In 2 (2) - 801 ln(2) + a 2 ) ■ 5, 



6 

5 

3( ai -41n(2))-^ + ^ 



r(n + l/2)r(n + 5/2)/ T/ _, T , 5. T . 

5 » = vf Z^r? ,1 - (*(» + 2) + *(n + l)-*(n+-)-*(n + 
1 (n + 2)1 (n + 1) V 2 
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where K = 0.915965 • ■ ■ is Catalan's constant and the other parameters, constants and 
series are: a x = 2/5-wi, a 2 = l-ir 2 -4iri/5, a 3 = -67r 2 /5 + 48193/7500 + 7r(7r 2 -3) i, 
i m = -0.740250494 ■ • •, r 53 = 2.22 5 2 4 6 651 • • and 

_ 1 - 6t + 2t 2 _ 3-12t + 8t 2 

2 2(t-l) ' 3 8(i- 1) 3 /2 

= [2, 41 /2, 313/48, 3047/480, ■ • •] , 
S% = [1, -25/2, -61/8, -129/16, ■ ••], 
5|° = [0,1, 7/10, 47/64, 981/1280, •••], 
SjjS = [0, 0, 161/300, 2039/4800, • • •], 
S"^ = [0, 0, 1847/18000, 2627/36000, • • •] , 
S$£ = [0, 0, 0, 14423879/7200000, • • •] 

The last bracket in H50|) comes from x^ 2 -*. 

Having the singular part of x^ a t the ferromagnetic and anti-ferromagnetic 
critical points, it is straightforward to obtain the asymptotic behavior of the series 
coefficients. This time, one needs the form of the coefficients in the expansion of 
ln 3 (l — x) that we find to be|| 

ln 3 (l-z)= £ ("-(*(»)+ 7) ( 51 ) 

n— 3 

where ^(l,n) is the first derivative of $>(n). Expanding ^(n) and ^(l,n) up to 1/n 2 
for large values of n, one obtains the following asymptotic behavior for the coefficients 
of the series: 

ln 2 (n) ln(n) 

c(n) ~ i 4 



128tt 3 n 128tt 3 n 2 

61 ln(n) 62 61 — 1 



64tt 3 n 2304tt 3 n 128tt 3 n 2 



where: 



h= 7 + 41n(2)-^, 
6 



b 2 = 288 In 2 (2) + 1447 ln(2) + I87 2 - 2IO7 - 8401n(2) + 45tt 2 +214 
7. x*- 1 -* + x^ versus x at scaling 

Thus far we have discussed, in Sections 4 and 6.2 the mathematical aspects of the 
solutions to the Fuchsian differential equations for x*- 3 -* and x^ ■ However, the physics 
implications of the solutions we have obtained call for some remarks near the physical 
critical points. Taking, as an example, the ferromagnetic singularity for x^ 3 \ the sum 
of the first two n-particle terms behave at r ~ as: 

m m l + lt ln 2 0) f , s 23\ ln(r) , . 

xW+x(3)~ —L- -_y + ln(2)-- --ii (52) 



16 7T 2 V 24 / 47r 

11 3, 1 / 2/ . 23, 14 

h -It o In (2) n(2 H 

48 8 3 4tt 2 V 12 1 ; 144 



|| An asymptotic form can be o btained using various packages available at 
http://algol .inria.fr/libraries/softwarc.html like the command "equivalent" in gfun 1241 . see 
details in 25 251. 
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The exact susceptibility, as reported in [16], yields for the normalized susceptibility 
X- 

-1/2 



- _ s (r+ \fi + t 2 )~ 

X ~ (l_ a 4)l/4 ' * - (l + r 2)l/8 X (53) 

— 1/4 oo oo 

(d r- 2 F + (r) + £ £ ■ - 9 ln P (r; 

where Ci = 1.000815260- • • is given with some 50 digits in JS]- -^+( T ) an d are 
given in ^3]. The constants 1 + l£ and ci verify 1 + 1^ + I§ = c\ with 9 digits, 7g~, 
corresponding to x^ 5 \ is the constant given in P (and with some 30 digits in 0). 
Thus, and as suggested in pp, the partial sums of the x would converge rapidly 
to the full x- Furthermore, adding x^ term has resulted in a series expansion that 
reproduces the first 24 terms of x to be compared with only eight first terms for x ' 
series. 

However, equation l|53|) shows a r -1 / 4 divergence as an overall factor to the 
logarithmic singularities. This structure, absent in Q52[). could suggest, in the most 
pessimistic scenario, that the n-particle sequence is perhaps useless in understanding 
scaling corrections and that one should be cautious in accepting the conclusions of 
studies of higher field derivatives of the susceptibility, based on similar n-particle 
representations fZ7\ |2E| • The same situation occurs for the low temperature regime 
when we compare the first two n-particle terms (x*" 2 "* and X ) with the full x & t 
scaling f . 

This observation raises several profound issues, which we do not address here. One 
is how the logarithmic terms in the entire sum add up to make the r -1 / 4 divergence be 
factored out. If one assumes that the other ^( 2n+1 ) terms share the same singularity 
structure as x 3 , m particular the occurrence (in variable r or s) of only integer critical 
exponents at the ferromagnetic critical point, the r -1 / 4 divergence, as an overall factor, 
implies the following correspondence : 

oo N(n) oo oo 

£ J2 • S n,m(r) ln m (r) - r^ 4 ■ £ £ b^ q) T« In"(r) 

n— 1 m— p—Oq—p 2 

with S n ^ m (r) analytical at r = and a n ,m numerical coefficients. N(n) is 
the maximum power of logarithmic terms occurring in the solution around the 
ferromagnetic point of the differential equation of x 2rl+1 ■ This correspondence requires 
probably a very particular structure in the successive differential equations. Obtaining 
the differential equation for x^ (or for x*- 6 -*), and obtaining much larger series for the 
full susceptibility Xj Wlu certainly help to guess such a structure and understand the 
susceptibility of the two-dimensional Ising model which continues to be a treasure- 
trove of profound insights into both the mathematics and physics of integrable systems. 

Let us note that the phenomenon we have discussed may be more widespread than 
that observed here. If so, a whole new chapter could be opened on field-theoretical 
expansions. The challenging problem one faces here is to link linear and non linear 
descriptions of a physical problem, namely the description in terms of an infinite 
number of holonomic (linear) expressions for a physical quantity of a non linear nature. 
Actually the latter is "Painleve like" since its series expansion can be obtained from 

% For the leading amplitude, x' 2 ) and x^' give 1/127T + 1^ ~ 1.0009593 ■ ■ ■ /12ir which is very close 
to 1.0009603 ■ ■ ■ /12n for the full % ED 
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a program of polynomial growth which uses exclusively a quadratic finite difference 
double recursion generalizing the Painleve equations |1 51 . The difficulty to link 
holonomic versus non-linear descriptions of physical problems is typically the kind 
of problems one faces with the Feynman diagram approach of particle physics, but 
the susceptibility of the Ising model is, obviously, the simplest non trivial example to 
address such an important issue. 

8. Conclusion 

We have introduced a simple and very efficient method to calculate numerically, 
with an arbitrary number of digits, the connection matrices between the independent 
solutions, defined at two singular points, of differential equations of quite high orders. 
We have considered the order seven, and ten, Fuchsian ODE's corresponding to the 
three and four particle contribution to the magnetic susceptibility of the Ising model. 
The entries of the connection matrix between two regular singular points have been 
obtained in floating point form and most of them have been recognized, particularly 
those that show up in the singular behavior of the physical solutions. They are 
expressed as polynomial, or algebraic, combinations of tt, ln(2), • • ■, radicals, and 
more involved numbers (not yet recognized) such as the " ferromagnetic constant" (J2J) . 
The method allows us to obtain the series expansions of the physical solutions x*- 3 -* 
(and x^) around any other regular singular point, besides the already known series 
around w = 0. We obtained, in this way, near each singular point all the dominant, 
and subdominant, singular behaviors of the physical solutions. Such subdominant 
singular behavior is certainly hard to obtain from series analysis. At the newly 
found quadratic singularities of the differential equation, we showed that the physical 
solution x' 3 ' itself is not singular. Also note, at w — 1/4, that the behavior in 
(1 — 4w)~ 3 / 2 corresponding to the largest critical exponent for the ODE is actually 
absent in the physical solution. Note the remarkable fact that the factorization of 
differential operator L7 (and £10) associated with x*- 3 ^ (respectively x^) shows clearly 
the differential operator responsible of the non-physical singularities given in E] 
and the newly found quadratic numbers 8 . In both cases (x an d x*- 4 '), these 
non-physical singularities are carried by the differential operator Z2 ■ N\ (respectively 
L25 • L12 ■ L 3 ■ L ) occurring at the right of L7 (respectively Cio). 

The physical solutions x' 3 -* (and x*- 4 )) being known as series around w = 0, the 
growth behavior of the corresponding series coefficients should be controlled by the 
singular behavior at the nearest singular points which are the ferromagnetic and anti- 
ferromagnetic critical points in both cases (10 = ±1/4 and x = 1). This growth 
is easily found from the expansion around the ferromagnetic and anti-ferromagnetic 
points. 

The connection matrices we have obtained allow us to relate the solutions around 
any given singular point to a common (non-local) basis of solutions. In this respect, 
we have obtained the exact expression of all the monodromy matrices, expressed in 
the same basis, and we have seen that they are simple matrices with rational function 
entries. In a forthcoming publication |14| . we will give the whole structure of the 
differential Galois group for the two previous Fuchsian differential equations. 

As far as the physics implications of the solutions are concerned, we have 
compared the corrections to scaling at the ferromagnetic point given by the first two 
terms (x^ and x^) with the full \. Qualitative difference is found raising profound 
issues on the n-particle representation of the susceptibility. The same observation 
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occurs for the antiferromagnetic point, and also for the low temperature regime. 
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9. Note added in the Proofs 

After completion of the revised version of our manuscript we were told that, as 
consequence of the work of B. M. McCoy, C. A. Tracy and T.T. Wu, the two 
transcendental numbers 1% and can actually be written in terms of poly logarithms, 
namely the Clausen function CI2 and of the Riemann zeta function, as follows : 

The derivation of these results has never been published but these results appeared 
in a conference proceedings [2H|- We have actually checked that J3 and we got 
from the calculations displayed in our paper as floating numbers with respectively 
42 f digits and 431 digits accuracy are actually in agreement with the previous two 
formula. These two results provide a clear answer to the question of how "complicated 
and transcendental" some of our constants occurring in the entries of the connection 
matrices can be. These extremely interesting results are not totally surprising when 
one recalls the deep link between zeta functions, polylogarithms and hypergeometric 
series [S01ET1E21I^1- 



10. Appendix A 

We give, in this Appendix, the explicit expressions of the differential operators X\ 
and Z2 and I3. The order one differential operator reads 

Ar=^+^ (54) 
aw pi 

with: 

pi = (-1 + w) (Aw - f) (1 + 2w) (Aw + 1) (l + 3w + 4w 2 ) 
(f - 3 w - 18 w 2 + 104 w 3 + 96 w 4 ) 

(l-7w- Aw 2 - A7w 3 + 36 w 4 + 280w 5 + 160w 6 + 256w 7 ) 
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p a = w- (-58 + 909 w + 3284 w 2 - 24711 w 3 - 72352 w 4 + 181016 w 5 

+ 1251768 w 6 + 2852880 w 7 + 1454592 w 8 - 11455616 w 9 

- 31712256 w 10 - 20418560 w 11 + 20840448 w 12 + 34963456 w 13 
+ 30146560 w 14 + 15728640 w 15 ) 

The order two differential operator Z 2 is 
1 2 rl 

Z 2 = - £ p» ■ A ( 55 ) 

where the polynomials Pi's, now, read: 

p 2 = w ■ (Aw - l) 2 (4 to + 1) (l + 3w + Aw 2 ) (-1 + w) (1 + 2w) 
(1 - 3w - I8w 2 + 104 w 3 + 96w 4 ) 

pi = (4w-l)(l-6u;-lllw 2 -108w 3 + 1080u; 4 - 4488w 5 

- 40368 w 6 - 94272 w 7 - 48384 w 8 + 72704 w 9 + 49152 w 10 ) 

p a = 4 + 48 w - 276 w 2 - 1520 w 3 - 3192 w 4 - 4224 - 71552 w e 

- 307200 w 7 - 239616 w 8 + 98304 w 9 + 98304 w 10 
The order three differential operator Y 3 is given by 

Y 3 = ±^ Pn .^ (56) 
P3 ^ dw" 

where the polynomials p^s, now, read: 

p 3 = w 2 ■ (w- l)(l + 2io) (l + 3w + Aw 2 ) (57) 

(Aw- if (Aw + if (96w 4 + 104 w 3 -18w 2 -3w+ if 

(l + 19 w - 368 w 2 - 3296 w 3 + 17882 + 272599 w 5 + 160900 w 6 

- 6979208 w 7 + 7550800 w 8 + 203094872 w 9 - 278920192 w 10 

- 3959814304 w 11 - 2115447424 w 12 + 20894729472 w 13 

+ 39719728128 w 14 + 20516098048 w 15 + 256763363328 w 16 

- 327065010176 w 17 - 8810227761152 w 18 + 414933057536 w 19 

+ 116411936538624 w 2a + 296827723186176 w 21 + 317648030138368 w 22 
+ 1 79148186189824 w 23 + 194933533179904 w 24 + 112931870081024 w 25 

- 55246164328448 w 26 + 1106 3 8 3 5 7 5 4 4 96 w 27 + 1511828488192 w 28 ) 

p 2 = w- (Aw - l) 2 (4w + l) 2 (96w 4 + 104w 3 - 18w 2 -3w + l) 2 
(6 + 102 w - 2018 w 2 - 23962 w 3 + 242904 w 4 + 2575633 w 5 

- 12389010w 6 - 178413527 w 7 + 80727412 w 8 + 6252221348 w 9 
+ 2456938016 w 10 - 178278888104 w 11 - 103902989696 w 12 

+ 3814815965856 w 13 + 1524977514176 w 14 - 67400886678400 w 15 

- 74115827788032 w 16 + 797710351468032 w 17 + 2324376661856256 w 18 
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- 1561280104050688 w 19 - 16314064973299712 w 20 

- 27005775986622464 w 21 - 40259640226480128 w 22 

+ 35764751009841152 w 23 + 1007304244270727168 w 24 
+ 1460771505523654656 w 2b - 13359756413056843776 w 26 

- 63988213537189134336 w 27 - 116684614339309600768 w 2S 

- 75710498024932245504 w 29 + 57121462326803824640 w 30 

+ 132479693600191414272 w 31 + 111232702128767107072 w 32 
+ 106152703871500156928 w 33 + 83508376521540632576 w 34 
+ 10084606300752183296 w 35 - 9404395631251816448 w 36 
+ 2682738003029262336 w 37 + 297237575406452736 w 38 J 

pi = 2- {Aw- l)(4io + l) (96u; 4 + 104w 3 - 18w 2 -3iu + l) 

^-3 - 25 w + 1013 w 2 + 7893 w 3 - 353904 w 4 - 1562671 w 5 

+ 43285825 w 6 + 192457911 w 7 - 2690351207 w 8 - 15077420736 w 9 
+ 94510776436 w 10 + 707838800508 w 11 - 2327528107216 w 12 

- 23421365465744 w 13 + 45755890012000 w 14 + 568028144875200 w 15 

- 824814656530816 w ie - 10390722028797440 w 17 

+ 12438134957505536 w 18 + 145637031330319360 w 19 

- 127616737495506944 w 20 - 1708173874007113728 w 21 

- 52355400373420032 w 22 + 15741676181476802560 w 23 

+ 24085046332129804288 w 24 - 57977682482294161408 w 25 

- 168033877030234750976 w 26 - 56941336876602621952 w 27 

- 426707803148891717632 w; 28 - 200805832817071095808 w 29 

+ 8716841486700848873472 w 30 - 6642009916749838811136 w 31 

- 192590979400145399971840 w 32 - 564260086660360537374720 w 33 

- 585770764250229243904000 u; 34 + 235172208485444226121728 w 35 

+ 1203159617695281059987456 w 36 + 1323272087085206269329408 w 37 
+ 997072075164663150542848 w 38 + 789138181323007857786880 w 39 
+ 388137877034203055390720 w 40 + 4946627729914186432512 w 41 

- 26947297377570617556992 w 42 + 10614515947351012540416 w 43 
+ 998718253365681192960 w; 44 ) 

pa = 2 w ■ (^-348 + 2768 w + 248784 w 2 - 358217 w 3 - 50461860 w 4 

+ 16394998 w 5 + 5283255372 w 6 + 3911764831 w 7 - 329364073508 w 8 

- 572985025996 w 9 + 13847002317264 w; 10 + 38091073842520 w 11 

- 437846238222272 w 12 - 1682624909395232 w 13 

+ 10892230218721408 w 14 + 52959188332189824 w 15 

- 214291413015639808 w 16 - 1200734422407578112 w 17 

+ 3319489124092462080 w 18 + 20066023020568346624 w 19 
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- 38248948302383529984 w 20 - 254480826931 185762304 w 21 

+ 261281404771497082880 w 22 + 2480194764802183397376 w 23 
+ 148352203759030894592 w 24 - 19049822668612433870848 w 25 

- 29328532357149024583680 w 26 + 103410036785394615320576 w 27 
+ 391034390334579595542528 u; 28 + 11096790708133489016832 w 29 

- 1530120948962096058466304 w 30 - 2868669407093825701150720 w 31 

- 6126661019209831555268608 w 32 + 2808943911875675603075072 w 33 

+ 40458568379798955017371648 w 34 - 169712327643359793079386112 w 35 

- 1092943871171162347998806016 w 36 - 1781375524629107822238367744 w 37 
+ 250471471742289487729786880 w 3S + 4679788548889591917580386304 w 39 
+ 7101176295364126941625974784 w 40 + 5918768536906007398653624320 u; 41 
+ 4083406571846803705271681024 w 42 + 2567747434748530216944009216 u; 43 
+ 846246487598480459424595968 w 44 - 49595159800068478383161344 w 45 

- 37040268890013610134208512 w 4e + 21784239691989525951676416 w 47 
+ 1753178556765355785584640 w 48 ) 

11. Appendix B: Solutions of the differential operator Y$ 

Considering the critical exponents at the regular singular points, as well as the formal 
solutions of differential operator Y3, one can make the following remarks. The roots 
of the polynomial of degree 28 in polynomial p^ (see H57[l) are apparent singularities. 
The roots of the polynomial of degree four in one of the factors of the same polynomial 
Ps are not apparent singularities. While the formal solutions near w = 0, w — ±1/4, 
and w — 00, have one Frobenius solution and two logarithmic solutions, the formal 
solutions near the other regular singular points are free of logarithmic solutions. The 
critical exponents at w = 1, w = —1/2, roots of 1 + 3w + Aw 2 — 0, and roots of 
1 - 3w~ 18 w 2 + 104 u; 3 + 96w 4 = 0, are respectively (-1,0,1), (-1,0,1), (-1,0,1) 
and (—1,1,2). This leads us to look for the solutions of the third order differential 
operator Y3 as a linear combination of powers of elliptic integrals with a common 
factor "taking care" of the non logarithmic singularity behavior of the singular points. 
Defining 

K(x) = 2 F 1 (1/2, 1/2; l;x) , E(x) = 2 F 1 (1/2,-1/2; l;x) 

and 

s(w) = w 2 ■ (l-16w 2 ) 3 (l + 2u;)(l-w)(l + 3w + 4u; 2 ) 
(l - 3w- 18w 2 + 104 w 3 + 96w 4 ) 
one obtains the three independent solutions of the differential operator Y 3 as: 

S!(Y 3 ) = J- . fp t . K 2 (l6w 2 ) + P 2 ■ E 2 (16w 2 ) 
s(w) V 

+ P 3 ■ K(l6w 2 )E{16w 2 ) 



S 2 (Y 3 ) = —— • P 4 • K\l/lQw z ) - \Qw z P 2 ■ E\l/lQw z ) 
s(w) 
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s(w) 



1 



+ P 5 • K(l/I6w 2 ) P(l/16w 2 )) 
y • ((Pi +P2 + Ps) ■ K 2 {1 - 16w 2 ) + P 2 
(2P 2 + P 3 ) • K(l - I6w 2 ) E(l - 16w 2 )) 



P 2 (l- 16w 2 ) 



with 



16 w 2 



P 



(1 - 16w 2 ) 2 
16 w 2 



1 - 16 
16 w 2 



P 5 = -2(1-16™ 2 )-P 2 -P 3 
where the three polynomials Pi, P 2 and P 3 read: 

Pi = - (l + 4w) (l - 5 u;-69 w 2 + 537 w 3 + 2964 w 4 -4100 w 5 

- 46816 w 6 - 74688 w 7 + 230656 w 8 + 647680 w 9 + 475136 w 10 
-8192w n + 720896 w 12 ) 

P 2 = - 1 + 5 to + 25 w 2 - 9 w 3 - 2408 w 4 - 17460 w 5 - 19696 w 6 

+ 28800 w 7 - 3328 w 8 - 62464 w 9 - 36864 w 10 
P 3 = 2 • (l - 3 w - 65 w 2 + 143 w 3 + 3888 w 4 + 15144 w 5 - 10624 w 6 

- 172416 w 7 - 241536 w 8 + 111616 w 9 + 282624 w w 
+ 180224 w n + 98304 w 12 ) 

Remark: Let us note the very close similarity between the differential operator 
Y3, occurring at the left of differential operator Lq (see Q) for x^ 3 \ and the differential 
operator M 2 (see JHJ) occurring at the left of differential operator £ 8 for . For 
this order four differential operator M 2 , we have been able, using the same ansatz, to 
obtain in closed form three of the four solutions, also expressed as a linear combination 
of products of elliptic integrals. Note that, setting A = 16w 2 , one can detect in the 
solutions of Y 3 (and also in the three solutions of M 2 we have found) the structure of 
S 3 permutation group |34) . A, 1/A, 1 — A, 1 — 1/A, etc. 

12. Appendix C: Connection matrices between w — and w = 1/4, w = 00 
12.1. Connection matrix between w — and w = —1/4 

The basis of solutions at the anti- ferromagnetic critical point w — — 1/4 are chosen as 
follows (with x = 1 + Aw) 




[1, 0, 1/10, -87/700, -313/1680, • • •], 
[0, 1, -17/10, -23/25, -1/30, • • •], 
[1, -5/2, -3/8, 5/16, 83/512, • • •], 
S { 4 ' 1/4) (x)-(Hx/8) + 6) + S^ 1/i) (x) 
Si~ m (x) • (ln 2 (x/8) + 121n(x/8) + 23264/315) 
+ 2Si- 1/4 \x) ■ (ln(x/8) + 6) + m {x) 
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with: 

S^ 1/4) (x) = [0,97/6,553/240, -2339/672, -1678457/645120, •••], 
s (k) 1/4) ( x ) = [0,0,0,85997/18000,8450503/1814400, •••]. 

Here again, an optimal choice of the components is made in order to remove logarithms 
and have as many zeroes as possible in the entries of the matrix. The same method of 
matching the series-solutions at a half-way point between w = and w = —1/4, gives 

C(0, -1/4) = (58) 



1 

















-2 


r-22 


r23 











-2iri 


r-32 + r 22 7r z 


^33 + 7*23 7T i 











6 


£*62 


£*53 








1 

8^ 


f + 6iri 


0,52 


053 





1 

16 






a 6 2 


063 


32 


-i 


1 

8 



with: 

r 22 r 3 3 ~ t-237-32 - 25/12288 



0152 = 


-3r 32 


5 . . 

- ^22 + «52«, 


«53 = ~3r 33 - 


5 . . 

^23 + *53« 


a-62 = 


,25 
' 16 ~ 


2tt 2 5tt n 


+ 67r r 32 


- 2-52 7T 


a-63 = 


,25 
' 16 ~ 


2tt 2 5tt n 

— - T 7)r 23 - 


+ 67r r 33 


- «53 7T 



and wheref: 

r 22 ~ -0.059050961331, r 23 ~ -0.018643190255, 

r 32 ~ 0.1631382423131, i 52 ~ -1.839621665835, 
i 53 ~ -0.015467563102 

12.2. Connection matrix between w = and w = oo 

The basis of solutions at the singular point w — oo are chosen as follows (with 
x = 1/w): 



s{°° 


(x) 






si°° 


(x) 


= [1,1, 7/16, 1/16, 7/256, •••], 




st 


(x) 


= (ln(x/4)-2/3).5(° o) (^) + ^ 


>(*)■ 


c(°° 


(x) 


= [0,1,0, 1/32, -9/512, •••], 




D 5 


(x) 


= (ln(.x/16) + a 1 )-5| oc) (x) + ^ 


(*). 


c(°° 
°6 


(x) 


= (ln 2 (x/16) + 2a 1 ln(:r/16) + a 2 ) 





+ 2(ln(x/16)+ ai ).^W + 5 6 rW 



| The numbers r;j's and j^'s are peculiar to each connection matrix. 
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with: 



Oi 



7T 2 379 r 
a 2 = — — + — + 57Ti, 



,30 w - [2/3, 1/6, 1/24, -1/96, 7/768, ■••], 

(x) = [0, 0, -3/2, 3/64, -107/512, -23113/491520, • • •] 
S^\x) = [0,0,0,93/44,-80891/13516,105811/4055040,- 
The connection matrix reads 
C(0, oo) = 



1 














1 


i 

16 


167T ' 








— 7T i 





1 

16 








11 + V4ii 


£42 - £i 


2 15 • 
tt 2 16tt 1 








051 


052 


9 49 ■ 

16 64tt 1 





1 

16 


aei 


062 


11 , 5tt • 75 • 
32 T 16 256tt 


64 


_ JLj 



where: 



^42 

a 5 i 



-1.534248223197, 

5 7T 

-J/41 + 77Ti, 



4 
29 
16 



y 41 ~ -22.932479960454, 

11 1 7T 

aB = -64-2 a!4a *-M*' 



+ ■ 



16^ 2 



7T 2 57T 

T ^4i + T 



062 



25 
256 



77^ 

192 



13. Appendix D 

13.1. Basis of solutions for w = 1, w = —1/2 and 1 + 3u> + 4w 2 = 0. 
The basis near w = 1 is (with x = 1 — w) : 

5(JVi)(a5), 

[0,0,0,1,65/24,383/72, •••], 
Sf ' (ln(x/24) + 2666/75) + S$(x), 
[0, 1, 0, 0, 0, -213149176769/914630737500, • 
[0, 0, 1, 0, 0, 806017240807/426827677500, • • 
[0, 0, 0, 0, 1, 555108887/158084325, • • •] , 



si 1 


(X 


si 1 


(X 


si 1 


(x 


si 1 


(x 


g(l 


(x 


c(l 
a 6 


(X 



with: 



,, tj ,., : - [0,96/5,628/25,0,-812657/18000,' 
The basis near w = —1/2 reads (with x = l + 2w) 



s[- 1/2 \x) 
St 1/2 \x) 



S(Wi)(aO, 

[0,0, 0,1, 8/3, 46/9, 247/27, •••], 
S(- 1/2) (x)-In(*) + ^ V2) (a;), 
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[0,1,0,0,0,-55489/60345, 




- 1/2) M = 


[0,0,1,0,0,159977/80460, 




- 1/2) M = 


[0,0,0,0,1,1492/447, •••] 



where: 



4o V2) (a) = [0, 3/4, 7/8, 0, -95/144, • • •] . 



The basis near w\ — — 3/8 + i y/7/8 root of 1 + 3w + Aw 2 is (with x = 1 — w/w\) 
S[ wi \x) = SiNiftx), 

Si wi) {x) = [0, 1, 49/64 - 61/(64 V7) i, 655/1024 - 747/(1024V7) i, ■ ■ ■ 



S^\x) 



: (m>i). 



«, -^r'O'-): ^ i} w= [0,0,1,0,0,. 

S$° l) (x)= [0,0, 0,0, !,•••] 



with: 



, 5 J (x) = [0,0,0,1,0,- ••], 



^o ^) = 657/896 + 61/(128V7)i, 41203/43008 + 1991/(6144\/7) v 



13.2. Connection matrices for w = 1, w = —1/2 and 1 + 3u> + 4w 2 = 
For the singular point w = 1, the connection matrix with w = reads 

C(0, 1) 



where: 



"A " 




A " 






where 




and 


B C 




B 




1 










4 + %2ii 
-2iri 



73 ■ 
144 1 

_7T\/3 

216 

r-42 + |r 52 i + t 



13 
2 



+ + »61« 



r 44 + z 44 i r 45 + i 45 i 

f *44 + »54« f «45 + *55* 



9*54 



9*55 



^46 + *46« 
f «46 + *56* 
9*56 



read 



73 

1447T 



773 

487T 



l*6i + (¥ + I)*2i+*5i* r 52 + |r 62 z-4M* ff§ 



ttVI,- _ 2575 „■ llTTV^ 
432 



573 , 7L- 
576tt T 18 



2573 
2304tt 



i 2 i ^ 1.838093775180, i 4 i ~ 4.136525226980, i 5 i - -8.13898927603 
i 61 ~ 20.74366088704, r 42 ~ 2.542631644752, r 52 ~ -0.01184208897 
r 62 ^ -4.87108777344, r 44 ~ 1.622875171987, r 45 ~ 1.954781507112 
r 46 ~ -3.51387499953, i 44 ~ 0.158271118920, i M ^ -2.13873967059 
i 45 ~ 0.041310289307, i 55 ^ -2.46759854730, i 46 ~ -0.02873064396 
i 56 ~ 4.392293882282, 



Galois group for the Ising model 
These numbers are such that: 

«46«55^44 + 745*56*44 ~ ^46*55*44 ~ *46 r 45'*54 

+ i45^46*54 - «45*56*'44 = - 



34 



468398 



18984375tt 2 

The connection matrix between w — and the singular point w = —1/2, reads 

A " 

where 



C(0, -1/2) 
1 



A 
B C 



B 



and [ C ] read 



_ 2 _ 31n(2) . 

-7T2 + 3 ln(2) 
r 4 i + uii 

f"51 + *51* 

r6i + «6i« 

f*44 + *44» 





^22 



ttV3 
54 

r-42 - 

»*52 - 



- 7T r 2 2 « 

- ? 42 i 

' *52* 



_V3 

97T 



ljV3 
9-rr 



5\/3 
36tt 



7"62 + 1621 

ta5 + iiz>i 



25yg 

1447T 



773 ■ 

9 4 

27 ~ f 



5V3 j 
18 



f46 + «46« 
3tt, 



-i 44 + i 54 i -^-145 + hzi -x*46 + *56* 



7"64 - — «44* r 65 - —l i5 l 



where: 



r 22 ~ -0.02539959775, r 4i ~ 6.805351589429, r 5 i ~ 7.203810787172, 
r 6 i ~ -8.75798651623, i 41 ~ -5.23529215352, i 5 i ~ 12.14972643902, 
z 6 i ^ 7.505979318469, r 42 ~ 0.512271205543, r 52 ^ -0.75497554989, 
r 62 = 2.232400538972, i 42 ^ 0.462196540081, i 52 ^ 0.143220115658, 
i 62 ~ -0.18195427623, r 44 ~ -0.1681290553, r 45 ~ -0.00270658055, 
r 46 ~ -0.00323043290, i 44 ~ -0.14301292413, i i5 ~ 0.690508507395, 
i 46 ~ -1.26354926677, i 54 ~ -0.34844554701, r 64 ~ 0.812327323812, 
i 55 ^ -0.50108648504, r 65 ^ 2.347957990666, i 5 6 ^ 1.132041888142, 
r 66 ~ -5.35056326640, 
The connection matrix between w = and the singular point w\ = —3/8 + i >/7 /8 



root of 1 + 3w + 4w 2 

C(0, wi) = 



0, reads 
A 





B C 



where 



A 
B 



and [ C ] read 



1 



T22 



2tt 



2tt 



»"31 + "T^Jl* + T« r 32 + ^-^22* 



275^7 „■ 
16384 ' 

623tt ■ 
24576 



r 4 i + 
rsi + h\i 
f6i + «6ii 



^42 + *42« 
7-52 + «52« 
T&l + »62« 



-§a 



(-! + 



'25 
v 16 



4 1 3 
7T 



2 *<)a 



Galois group for the Ising model 



35 



where: 



r 44 + z 44 i r 45 + i 45 i r 46 + i 46 i 

7*54 + «54« ?"55 + «55« f 56 + «56* 
7*64 + «64« 7*65 + «65« 7*66 + «66« 



_ 825V7- 1869i 
° ~~ 16384tt ' 
r 2 i ~ -0.30983963151, r 3 i 

r 32 ~ 0.044743829620, r 41 : 
r 42 ~ 0.028522637766, z 42 c 
i 5 i ~ 10.77185269595, r 52 = 
r 6i ~ -6.98898250954, i 61 
i 62 ~ 1.339914984659, r 44 c 
r 45 ~ -0.02716280332, i 45 
i 46 ~ 0.050869745772, r 54 c 
r 55 ~ 0.496361798471, i 55 z 
i 56 ~ 0.040697038977, r 64 c 
r 65 ~ -0.14739127007, i 65 
i 66 ~ -1.29483325656, 



i 1.38629436111, r 22 : 

: 4.70316610599, i 41 - 
0.03731267544, r 5 i ^ 
0.25654299002, i 5 2 = 

t -17.585497074, r 62 
0.00394832042, i 44 - 

t -0.0900753899, r 46 
-0.2122947699, i 54 i 
0.00455966493, r 56 - 
-0.1279407612, z 64 c 

i 1.64596123266, r 66 : 



i -0.07996746793, 
-5.10203220992, 
1.404170417754, 
-0.03695328252, 

- -0.18342705750, 
0.043931830095, 

- 0.070134204478, 
: 0.033562029788, 
-0.36867647137, 

: -0.68382860060, 
; 0.189914623980, 



14. Appendix E: Monodromy matrices in the w = 0-basis 

The monodromy matrix around w = expressed in terms of its own (w = 0) well- 
suited basis is given in J33J). 

The monodromy matrix around w — —1/2, expressed in terms of the (w — 0) 
well-suited basis, after a conjugation similar to l|34l) . and thus using the previously 
given connection matrices, reads in terms of a and ft: 

A 



where: 



4a 2 ■ M w = (-l/2){a, ft) = 



4 a 2 
48 a ft 



B C 



A 
B 



24 a 2 
-528 a ft 
-12 (14 a + 5) aft 
—a a ft 





4a (12ft- 
24 ft a 2 
-528aft 



•a 





-96 ft 
4 (a - 12 ft) a 
1056 ft 



-12 (14a + 5) aft 24(14a + 5)ft 
—a aft 2 ft a 



with a = (-75 + 52 a 2 + 60 a) and [ C ] = 4a 2 • Id(3 x 3). 

The monodromy matrix around w = 1/4, expressed in terms of the (w — 0)-well 
suited basis reads: 

" A 

24a 4 - M„ =0 (l/4)(a, ft) = 

B C 
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A 



36 



where 



B 



and [ C ] read respectively: 



-24 a 4 
-48 a 4 






24 a 4 






-144 a 2 
24 a 4 



-48 (5 a 4 + 8 O 2 + 8 2 a 2 ) 32 (4 a 2 - 75 - 15 a 2 ) O 48 (9 a 2 + 80 O) O 
12 (5a 2 + 40 + 40a 2 ) a 2 4 (75 - 4 a 2 ) a 2 -300a 2 O 

-(87 + 8a 2 )a 4 3(4a 2 -75)a 2 



and: 



24 a 4 





-384 a 2 1536 O 2 
24 a 4 -192 a 2 
24 a 4 



The monodromy matrix around w = —1/4, expressed in terms of the (w = 0) 
well-suited basis reads: 



where 



A 
B 



12a 4 -M„ =0 (-l/4)(a, O) 

and [ C ] read respectively: 




A 
B C 



-12 a 4 
48 a 4 
24 a 5 
a 4 i 

051 

a 6 i 



12 a 4 


a 4 2 

0,52 

a-62 






12 a 4 
192 ^ (lOO-3a 2 ) 
48aO (-5a + 20O -6a 2 ) 
48a 2 (-5a+10O-3a 2 ) 



with: 



a 41 = -144 a 4 - 192 ft 2 - 192 2 a 2 , 

a i2 = -16 (60aO + 75O + 8Oa 2 - 18 a 3 + 15 a 2 ) , 

a 5 i = -12 a (5 a 3 + 6 a 4 + 8 O 2 + 8 2 a 2 - 2 a - 2 a 3 0) , 

a 52 = -2 a (300 + 32 a 2 + 240 a - 80 a 3 - 75 a) , 

a 6 i = -a 2 (-69 a 2 + 60 a 3 + 34 a 4 + 48 O 2 + 48 2 a 2 - 24 a O - 24 a 3 0) 



a 62 = -2a 2 (150 O- 30 a 2 



16 0a 2 + 120 aO -44a 3 - 75 < 



and: 



12 (q + 40) 2 a 2 



-192 (a + 40) aO 



768 O 2 



24a 3 0(a + 4 0) -12 a 2 (-a 2 + 32 O 2 ) 96 (-a + 4 0) a 
48 a 4 2 48 (a - 4 0) a 3 12 (-a + 4 O) 2 a 2 
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The monodromy matrix around w 
suited basis reads: 



oo, expressed in terms of the w 



37 
O-well 



where 



A 
B 



with: 



and: 



24a 4 • M w=0 (oo)(a, O) = 
and [ C ] read respectively: 



A 
B 





c 



24 a 4 








-288 a 3 


-24 a 3 (-a + 60) 


-864 a 2 


48 a 4 


24 a 4 


24 a 3 (a + 6 0) 


an 


0,42 


96 (-21a 2 + 160 0) 


asi 


«52 


-120 (-6 a 2 - a + 32 0) aO 


061 


«62 


6 (20 a - 225 - 36 a 2 + 160 0) a 2 



a 41 = 96 (a 3 + 16 a 2 - 16 O) O 

«42 = 16 (-33 a 3 - 60 a 2 + 240 a O + 8 O a 2 - 600 O) O, 

a 51 = -24 (2 a 3 - 15 a 2 - 4 a + 16 O a 2 - 16 0) a 0, 

a 52 = -4 (-40 a 3 -45a 2 - 150 a + 240 aO + 8 a 2 - 600 O) aO, 

a 61 = 6 (-20 a 2 - 83 a + 4 a 3 + 16 O a 2 - 16 O) a 2 0, 

a 62 = (-525 a - 44 a 3 + 240 a O + 8 O a 2 - 600 O) a 2 0, 



24(-a + 4 0) 2 a 2 768 (-a + 4 0) a 



6144 O 2 



-24a 3 (-a + 4 0) -24 a 2 (-a 2 + 32 O 2 ) -384aO(a + 4 0) 



24a 4 2 



!a 3 (a + 4 0) 



24 a 2 (a + 4 0) 2 
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